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Abstract 

We study the critical properties of the weakly disordered p-component ferro- 
magnet in terms of the renormalization group (RG) theory generalized to take into 
account the replica symmetry breaking (RSB) effects coming from the multiple 
local minima solutions of the mean- field equations. It is shown that for p < 4 the 
traditional RG flows at dimensions D = 4 — e, which are usually considered as de- 
scribing the disorder-induced universal critical behavior, are unstable with respect 
to the RSB potentials as found in spin glasses. It is demonstrated that for a general 
type of the Parisi RSB structures there exists no stable fixed points, and the RG 
flows lead to the strong coupling regime at the finite scale ~ exp(l/u), where 
u is the small parameter describing the disorder. The physical concequences of 
the obtained RG solutions are discussed. In particular, we argue, that discovered 
RSB strong coupling phenomena indicate on the onset of a new spin glass type 
critical behaviour in the temperature interval r < t* ~ exp(— — ) near T c . Possible 
relevance of the considered RSB effects for the Griffith phase is also discussed. 
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1 Introduction 



In this paper we study the effects produced by weak quenched disorder on the critical 
phenomena in the ferromagnetic spin systems near the phase transition point. In the 
most general terms the traditional point of view on this problem could be summarized 
as follows. 

According to the usual scaling theory near the critical temperature T c the only 
relevant scale that remains in the system is the correlation length R c which scales as 
~ t~ u , where r = (T — T c )/T c << 1 is the reduced temperature parameter and v is the 
correlation length critical exponent. 

If the disorder is weak (e.g. the concentration of impurities is small), its effect on 
the critical behavior in the vicinity of the phase transition point T c remains negligible 
so long as the correlation length R c is not too large, i.e. for temperatures T not too 
close to T c . In this regime the critical behavior will be essentially the same as in the 
pure system. 

However, in the close vicinity of the critical point, at r = (T — T c )/T c — > 0, the 
correlation length R c grows and becomes larger than the average distance between the 
impurities, so that the effective concentration of impurities, measured with respect to 
the correlation length, becomes large. The strength of disorder, described by small 
parameter u, affects only the width of the temperature region near T c in which the 
effective concentration gets large. If uR® ^> 1, where D is the spatial dimensionality, 
one has no grounds, in general, for believing that the effect of impurities will be small. 

A very simple general criterion has been discovered, the so-called Harris criterion 
which makes it possible to predict the effect of impurities qualitatively from only the 
critical exponents of the pure system. According to this criterion the impurities change 
the critical behavior only if a, the specific heat exponent of the pure system, is greater 
than zero (i.e. the specific heat of the pure system is divergent at the critical point). 
According to the traditional point of view, when this criterion is satisfied, the disorder 
becomes relevant and a new universal critical behavior, with new critical exponents, is 
established sufficiently close to the phase transition point p], |3|: 



This argument identifies l/« as the cross-over exponent associated with randomness 
[§]. In contrast, when a < (the specific heat is finite), the disorder appears to be 
irrelevant, i.e. their presence does not affect the critical behavior. 

Near the phase transition point the D- dimensional Ising-like systems can be de- 
scribed in terms of the scalar field Ginsburg-Landau Hamiltonian with a double-well 
potential: 



Here the quenched disorder is described by random fluctuations of the effective tran- 
sition temperature 5r(x) whose probability distribution is taken to be symmetric and 
Gaussian: 



T <T U =U 



1/a 
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P[5r] = Po exp(-±-Jd D x(6T(x)) 2 j , (1.3) 

where u <C 1 is the small parameter which describes the disorder, and po is the normal- 
ization constant. In Eq. ( |1 . 2| ) r ~ (T — T c ) and for notational simplicity, we define the 
sign of 5t(x) so that positive fluctuations lead to locally ordered regions. 

Now, if one is interested in the critical properties of the system, one has to integrate 
over all local field configurations up to the scale of the correlation length. This type 
of calculation is usually performed using a Renormalization Group (RG) scheme, which 
self-consistently takes into account all the fluctuations of the field on length scales up 
to R c . 

To derive the traditional results for the critical properties of this system discussed 
above one can use the usual RG procedure developed for dimensions D = 4 — e, where 
e C 1. Then one finds that in the presence of the quenched disorder the pure system 
fixed point becomes unstable, and the RG rescaling trajectories are arriving to another 
(universal) fixed point g* ^ 0; ^ 0, which yields the new critical exponents describing 
the critical properties of the system with disorder. 

However, there exists an important point which missing in the traditional approach. 
Consider the ground state properties of the system described by the Hamiltonian ( |1.2| ). 
Configurations of the fields 0(x) which correspond to local minima in H satisfy the 
saddle-point equation: 

- A(f)(x) + (r - 6r(x))<j)(x) + g<f(x) = . (1.4) 

Clearly, the solutions of this equations depend on a particular configuration of the 
function 5t(x) being inhomogeneous. The localized solutions with non-zero value of 
exist in regions of space where r — 5t(x) has negative values. Moreover, one finds 
a macroscopic number of local minimum solutions of the saddle-point equation ( |1.4| ). 
Indeed, for a given realization of the random function 5t(x) there exists a macroscopic 
number of spatial "islands" where r — 5t(x) is negative (so that the local effective 
temperature is below T c ), and in each of these "islands" one finds two local minimum 
configurations of the field: one which is "up", and another which is "down". These local 
minimal energy configurations are separated by finite energy barriers, whose heights 
become larger as the size of the "islands" are increased. 

The problem is that the traditional RG approach is only a perturbative theory in 
which one treats the deviations of the field around the ground state configuration, and 
it can not take into account other local minimum configurations which are "beyond 
barriers" . This problem does not arise in the pure systems, where the solution of the 
saddle-point equation is unique. However, in a situation like that discussed above, when 
one gets numerous local minimum configurations separated by finite barriers, the direct 
application of the traditional RG scheme may be questioned. 

In a systematic approach one would like to integrate in an RG way over fluctuations 
around the local minima configurations. Furthermore, one also has to sum over all these 
local minima up to the scale of the correlation length. In view of the fact that the 
local minima configurations are defined by the random quenched function Sr(x) in an 
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essentially non-local way, the possibility to implement successfully such a systematic 
approach seems rather hopeless. 

On the other hand there exists another technique which has been developed specif- 
ically for dealing with systems which exhibit numerous local minima states. It is the 
Parisi Replica Symmetry Breaking (RSB) scheme which has proved to be crucial in the 
mean-field theory of spin-glasses (see e.g. 0). Recent studies show that in certain cases 
the RSB approach can also be generalized for situations where one has to deal with 
fluctuations as well 0,0, [|J. 

It can be argued that the summation over multiple local minima configurations in 
the present problem could provide additional non-trivial RSB interaction potentials for 
the fluctuating fields |§. Let us consider this point in some more details. 

To carry out the appropriate average over quenched disorder we can use the standard 
replica approach. To do this, we need to average the nth (n — > 0) power of the partition 
function. This is accomplished by introducing the replicated partition function, Z n = 
Z u [5t], where (...) denotes the averaging over St{x) with the probability distribution 
( |1.3| ). Simple integration yields: 



Z n = Z^{8r) = JD(j> a {x)exp 



- / d D x | £™=i [V0 a (x)] 2 + \t E: =1 ft{x) 



1.5) 



where 



9ab = g5 a b ~ u 



1.6) 



is the replica symmetric (RS) interaction parameter. If one would start the usual RG 
procedure for the above replica Hamiltonian (as it is done in the traditional approach), 
then it would correspond to the perturbation theory around the homogeneous ground 
state = 0. 

However, in the situation when there exist numerous local minima solutions of the 
saddle-point equation (|1.4|) one has to be more careful. Let us denote the local solutions 
of the eq. (|1.4j) by ip^\x) where i = 1,2,... iVo labels the "islands" where 8r(x) > r. If 
the size Lq of an "island" where (Sr(x) — r) > is not too small, then the value of i/j^'(x) 
in this "island" should be ~ ±yj (5t(x) — r)/g, where St(x) should now be interpreted 
as the value of St averaged over the region of size L Q . Such "islands" occur at a certain 
finite density per unit volume. Thus the value of N is macroscopic: N Q = nV, where 
V is the volume of the system and k is a constant. An approximate global extremal 
solution $(x) is constructed as the union of all these local solutions without regard for 
interactions between "islands." Each local solution can occur with either sign, since we 
are dealing with the disordered phase: 



[x; St(x)] 



nV 
i=l 



X) 



;i-7) 
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where each <jj = ±1. Accordingly, the total number of global solutions must be 2 kV . We 
denote these solutions by $( a )[ar; 8r(x)), where a = 1, 2, K = 2 kV . As we mentioned, 
it seems unlikely that an integration over fluctuations around <fi(x) = will include 
the contributions from the configurations of <f)(x) which are near a since Q(x) is 

"beyond a barrier," so to speak. Therefore, it seems appropriate to include separately 
the contributions from small fluctuations about each of the many $( a Jx; 5r\. Thus we 
have to sum over the K global minimum solutions (non-perturbative degrees of freedom) 
St] and also to integrate over "smooth" fluctuations (p(x) around them 



Z[5r] = J Dip(x) J2a ex P ~H 



T 



(1.8) 



/ D(p(x) exp —H[ip; St] x Z[ip; 5t] 



where 



Z[<p;6T]=TZe X p[-H a -fd D 



x 



g$ 2 (a)V 2 (x)+g<5> {a)V 3 (x) 



1.9) 



and H a is the energy of the a-th solution. 

Next we carry out the appropriate average over quenched disorder, and for the replica 
partition function, Z n , we get: 

Z n = J D5r J D^exp J d D x[5r(x)} 2 - fj H[ Va , 6r}^ x Z n [ Va ,5r} , (1.10) 

where the subscript a is a replica index and 



z n [t£ a ; 5t] = ex p (- H H <*a- J d ° x 



9®\ aa ){x)V 2 a {x) + 9®(a a )(x)<pl(x) 



(1.11) 

It is clear that if the saddle-point solution is unique, then from the eq. (|1.10|) , (|l.ll|) 
one would obtain the usual RS representaion ( |1.5|) , (|L~6f) . However, in the case of the 
macroscopic number of the local minima solutions the problem is getting highly non- 
trivial. This sutuation is reminiscent of the (unsolved) problem of summing over the 
saddle-point solutions in the random-field Ising model, which is believed to provide the 
RSB phase near the phase transition point 
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It is obviously hopeless to try to make a systematic evaluation of the above replicated 
partition function. The global solutions are complicated implicit functions of 5t(x). 
These quantities have fluctuations of two different types. In the first instance, they 
depend on the stochastic variables 5r(x). But even when the 5r(x) are completely 
fixed, $( Q )(x) will depend on a (which labels the possible ways of constructing the 
global minimum out of the choices for the signs {a} of the local minima). A crude way 
of treating this situation is to regard the local solutions ip^(x) as if they were random 
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variables, even though St(x) has been specified. This randomness, which one can see is 
not all that different from that which exists in a spin glass, is the crucial one. It can be 
shown then, that due to the interaction of the fluctuating fields with the local minima 
configurations (the term <&^ a ^(p\ in the eq.( |l.ll|) ), the summation over solutions in the 

replica partition function Z n [(p a ], eq.( |l.ll| ), could provide the additional non-trivial RSB 
potential J2ab dab^Pa^l m which the matrix g ab has the Parisi RSB structure [|9] . 

In the paper |§ due to several simplifying assumptions, the matrix g a b has been 
obtained to have explicit 1-step RSB structure, which, in general, may not be the case. 
In this paper we are going to study the critical properties of weakly disordered systems 
in terms of the RG approach taking into account the possibility of a general type of the 
RSB potentials for the fluctuating fields. The idea is that hopefully, like in spin-glasses, 
this type of generalized RG scheme self-consistently takes into account relevant degrees 
of freedom coming from the numerous local minima. In particular, the instability of 
the traditional Replica Symmetric (RS) fixed points with respect to RSB indicates that 
the multiplicity of the local minima can be relevant for the critical properties in the 
fluctuation region. 

It will be shown in the next Section that, whenever the disorder appears to be 
relevant for the critical behavior, the usual RS fixed points (which used to be considered 
as providing new universal disorder-induced critical exponents) are unstable with respect 
to "turning on" an RSB potential. Moreover, it will be shown that in the presence of 
a general type of the RSB potentials the RG flows actually lead to the strong coupling 
regime at the finite spatial scale R* ~ exp(l/-u) (which corresponds to the temperature 
scale ~ exp(— -)). At this scale the renormalized matrix g ab develops strong RSB, 
and the values of the interaction parameters are getting non-small. 

Usually the strong coupling situation indicates that sertain essentially non-perturbative 
exitations have to be taken into account, and it could be argued that in the present model 
these are due to exponentially rare "instantons" in the spatial regions, where the value 
of St(x) ~ 1, and the local value of the field (f(x) must be ~ ±1. (Distant analog of 
this situation exists in the 2D Heisenberg model where the Poliakov renormalization 
develops into the stong coupling regime at a finite (exponentially large) scale which is 
known to be due to the non- linear localized instanton solutions |TT|). 

In Section 3 we discuss the physical concequences of the obtained RG solutions. In 
partiqular we show that due to the absence of fixed points at the disorder dominated 
scales R » u~ u / a (or at the corresponding temperature scales r « u l l a ) there must 
be no simple scaling of the correlation functions or of other physical quantities. Besides, 
we demonstrate, that the structure of the SG type two-points correlation functions is 
characterized by the strong RSB, indicating on the onset of a new type of the critical 
behaviour of the SG nature. 

In Section 4 we consider the special case of systems with the number of spin com- 
ponents p = 4, in which the pure system specific heat critical exponent a = 0. Here 
the disorder appears to be marginally irrelevant in a sense that it does not change the 
critical exponents. Nevertheless, the critical behaviour itself (described in terms of the 
logarithmic singularities) is effected by the disorder, and moreover, the RSB phenomena 
is demonstrated to be relevant in this case as well. 
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The remaining problems as well as future perspectives are discussed in the Conclu- 
sions. In particular, we discuss possible relevance of the considered RSB phenomena for 
the Griffith phase which is known to exist in a finite temperature interval near T c [[TIJ . 



2 Replica Symmetry Breaking in the Renormaliza- 
tion Group Theory 

We consider the p-component ferromagnet with quenched random effective temperature 
fluctuations, which near the transition point can be described by the usual Ginzburg- 
Landau Hamiltonian: 



H[5r,<f>] = Jd D x 



(2.1) 



+1(t - 5t(x)) ELi + \g Ef )i= i 



where the quenched random temperature Sr(x) is described by the Gaussian distribution 

(ED- 

In terms of the standard replica approach after integration over 8r{x) for the replica 
partition function one gets: 



fD$(x)exp 



fd D x UE?=iE:=i[V^(a;)] 2 + |rELiE^i[^(x)] s 



;2.2) 



where 



9ab = g^ab - u . 



(2.3) 



To study the critical properties of this system we use the standard RG procedure 
developed for dimensions D = 4 — e, where e<l. Along the lines of the usual rescaling 
scheme (see e.g. |L3||) one gets the following (one-loop) RG equations for the interaction 
parameters g a b- 



dg, 



ab 



di 



tg a b - -z-^(4:gl b + 2(g aa + g bb )g ab + p ^ g ac gcb) 
87r c =i 



(2.4) 



where £ is the standard rescaling parameter. 

Changing g ab — > 8n 2 g ab , and g a ^ b — > —g a ^ b (so that the off-diagonal elements would 
be positively defined), and introducing g = g aa , we get the following RG equations: 
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dg t n 

-7T = £9ab - (4 + 2p)gg ab + Ag 2 ab + p 9ac9cb (a^b), (2.5) 

c+a,b 

j n 

~g = e~g-(8 + p)g 2 - P J2gl (2.6) 

If one takes the matrix g a b to be replica symmetric, as in the starting form of Eq. 
( |2.3| ), then one would recover the usual RG equations for the parameters g and u, and 
eventually one would obtain the well known results for the fixed points and the critical 
exponents JJ, |3| . Here we leave apart the question as to how perturbations out of the RS 
subspace could arise (see discussion in ||) and formally consider the RG eqs.( [2.5|) , ([2.6|) 
assuming that the matrix g a b has a general Parisi RSB structure. 

According to the standard technique of the Parisi RSB algebra (see e.g. H), in the 
limit n — > the matrix g ab is parametrized in terms of its diagonal elements g and the 
off-diagonal function g(x) defined in the interval < x < 1. All the operations with the 
matrices in this algebra can be performed according to the following simple rules (see 
e.g. 1,0): 

<&-(sr*; (2-7) 

n 

(g 2 )ab = J2 9ac9cb -> (c; c(x)), (2.8) 



where 



c=l 

c = 9 2 - Jodxg 2 (x), 



c(x) = 2(g - Jq dyg(y))g(x) - ^ dy[g(x) - g(y)} 2 . 

The RS situation corresponds to the case g(x) = const independent of x. 
Using the above rules from the eqs.([2.5|), (|2.6| ) one gets: 



;2.9) 



—g{x) = {e-{A + 2p)~g)g{x)+Ag 2 {x)-2 P g{x) f dyg{y)-p f ' dy(g(x) - g(y)) 2 (2.10) 
at; Jo Jo 

-^~g = e~g-(8+p)~g 2 +p7 (2.11) 

where g 2 = Jq 1 dxg 2 (x). 

Usually in the studies of the critical behaviour one is looking for the stable fixed- 
points solutions of the RG equations. The fixed-point values of the of the renormalized 
interaction parameters are believed to describe the asymptotic structure of the effective 
Hamiltonian which makes possible to calculate the singular part of the free energy, as 
well as the other thermodynamic quantities. 

From the eq.( [2.1U ) one can easily find out what should be the structure of the 



function g(x) at the fixed point, -^g(x) = 0, -^g = 0. Taking the derivative over x 
twice, one gets, from Eq. ( |2.10|) : g'(x) = 0. This means that either the function g(x) is 
constant (which is the RS situation), or it has the step-like structure. It is interesting to 
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note that the structure of fixed-point equations is similar to those for the Parisi function 
q(x) near T c in the Potts spin-glasses |15|], and it is the term g 2 (x) in Eq. fl2.10p which 



is known to produce lstep RSB solution there. The numerical solution of the above 
RG equations convincingly demonstrates that whenever the triel function g(x) has the 
many-step RSB structure, it quickly developes into the 1-step one with the coordinate 
of the step being the most right step of the original many-step function. 
Let us consider the 1-step RSB ansatz for the function g(x): 

q(x) = i 9o for < x < x Q 
v ' \ gi for x < x < 1 v ; 

where < Xq < 1 is the coordinate of the step. 

In terms of this ansatz the above fixed-point equations have several non-trivial solu- 
tions: 

1) The RS fixed-point which corresponds to the pure system: 

go = gi = 0; g = -^—t (2.13) 
8 + p 

This fixed point (in accordance with the Harris criterion) is stable for the number of 
spin components p > 4, and it is getting unstable for p < 4. 

2) The disorder- induced RS fixed point (for p > 1)0, ||: 

A ~P ~ P (oia\ 

9o = 9l = e i6(^Ty 9 = e m(^Ty (214) 

It was usually considered to be the one which describes the new universal critical be- 
haviour in systems with impurities. This fixed point has been shown to be stable (with 
respect to the RS deviations!) for p < 4, which is consistent with the Harris criterion. 
(For p = 1 this fixed point involves an expansion in powers of (e) 1//2 and this structure 
is only revealed within a two-loop approximation). However, the stability analysis with 
respect to the RSB deviations shows that this fixed point is always unstable p. There- 
fore, whenever the disorder is relevant for the critical behaviour, the RSB perturbations 
must be getting the dominant factor in the asymptotic large scale limit. 

3) The 1-step RSB fixed point 0: 

#0 = 0; gi = e 
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16(p-l)-px (8+p)' 

(2.15) 



16(p-l)-px (8+p) ' 

This fixed point can be shown to be stable (within 1-step RSB subspace!) for: 

1 < p < 4, 

(2.16) 

< z < x c {p) = ifgf 
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In particular, x c {p = 2) = 4/5; x c {p = 3) = 32/33, and x c (p = 4) = 1. Using the 
result fl2.15|) one can easily obtain the corresponding critical exponents which are now 
getting to be non-universal being dependent on the starting parameter xq || (see also 
next Section). 

(Note, that in addition to the fixed points listed above there exist several other lstep 
RSB solutions which are either unstable or unphysical.) 

The problem, however, is that if the parameter xq of the starting function g(x; £ = 0) 
(or, more generaly, the coordinate of the most right step of the many-steps starting 
function) is beyond the stability interval, such that x c {p) < x < 1, then there exist no 
stable fixed points of the RG eqs.( [2.10| ),( [2.11|) . One faces the same situation, of course, 



in the case of a general continuous starting function g(x; £ = 0). Moreover, according to 
eq. (|2.16| ) there exist no stable fixed points out of the RS subspace in the most intersting 



Ising case, p = 1. 

Unlike the RS situation for p = 1, where one finds the stable ~ ^fk fixed point 
in the two-loop RG equations here adding next order terms in the RG equations 
doesn't cure the problem. In the considered RSB case one finds that in the two-loops 
RG equations the values of the parameters in the fixed point are formally getting of the 
order of one, and it signals that we are entering the strong coupling regime where all 
the orders of the RG are getting relevant. 

Nevertheless, to get at least some information about the physics behind this insta- 
bility phenomena, one can proceed analizing the actual evolution of the above one-loop 
RG equations. The scale evolution of the parametrs of the Hamiltonian would still ade- 
quately describe the properties of the system until we reach a critical scale £*, at which 
the strong coupling regime begins. 

The evolution of the renormalized function g(x; £) can be analyzed both numerically 
and analytically. It can be shown (see Appendix A) that in the case p < 4 for a general 
continuous starting function g(x; £ = 0) = go{x) the renormalized function g(x; £) tends 
to zero everywhere in the interval < x < (1 — A(£)), while in the narrow (scale 
dependent) interval A(£) near x = 1 the values of of the function g(x; £) grow: 

f a^; at (1 - x) « A(£) 

#(z;£)~ (2.17) 

[ 0; at (1 -x) » A(£) 

$(£)~uln— !— (2.18) 
1 — it£ 

where 

A(£)~(l-u£) (2.19) 

Here a is a positive non- universal constant, and the critical scale £* is defined by the 
condition that the values of the renormalized parameters are getting of the order of one: 
(1 — w£*) ~ u, or £* ~ 1/u. Correspondingly, the spatial scale at which the system is 
entering the strong coupling regime is: 

R, ~ exp(^) (2.20) 
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Note that the value of this scale is much bigger than the usual crossover scale ~ u~ a l v 
(where a and v are the pure system specific heat and the correlation length critical 
exponents), at which the disorder is getting relevant for the critical behaviour. 

According to the above result, the value of the narrow band near x = 1 where the 
function <?(x; £) is formally getting divergent is A(£) ~ (1 — u£) — > u « 1 as £ — > £*. 

Besides, it can also be shown (Appendix A) that the value of the integral g(£) = 
Jo 9( x i0 is formally getting divergent logarithmically as £ — > £*: 

^(0~Mln-^- (2.21) 
1 - < 

Qualitatively similar asymptotic behaviour for g(x;l;) is obtained for the case when 
the starting function poO^) nas the 1-step RSB structure (|2.12| ), and the coordinate of 
the step x Q is in the instability region (or for any x in the Ising case p — 1): 



gi(0) 



: ^ clt Xq <C X <C 1 



l-(4-2p+pxo)<?i(0)C' 

9&Z)~{ (2-22) 
0; at < x < xq 

Here (?i (0) = gi(£ = 0) ~ w, and the coefficient (4 — 2p + pxo) is always positive. In this 
case again, the system arrives into the strong coupling regime at scales £ ~ 1/u. 

Note that the above asymptotics do not explicitely involve e. Actually, the role of the 
parameter e > is to "push" the RG trajectories out of the trivial Gaussian fixed point 
g = 0; g = 0. Thus, the value of e, as well as the values of the starting parameters go(x), 
go, define a scale at which the solutions finally arrive to the above asymptotic regime. 
In the case e < (above dimensions 4) the Gaussian fixed point is stable; on the other 
hand, the strong coupling asymptotics still exists in this case as well, separated from 
the trivial one by a finite (depending on the value of e) barrier. Therefore, although 
infinitely small disorder remains irrelevant for the critical behaviour above dimensions 
4, if the disorder is strong enough (bigger than sertain depending on e threshould value) 
the RG trajectories could arrive to the above strong coupling regime again. 

3 Scaling and Correlation functions 

3.1 Temperature Scales 

The renormalization of the mass term T(£)^a=i^a i s described by the following RG 
equation: 

| lnr = 2 -^K 2 +^ + ?'E 9la) (3.1) 

Changing (as in the previous section) g ab — > 8n 2 g ab , and g a ^ — > —g a ^b, in the Parisi 
representation we get: 

-| In r = 2 - [(2 + p)g(£) + p £ g(x; £)] (3.2) 
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or 



r(£)=T exp{2£- f dr)[{2 + p)g{r)) + pg{r]))} (3.3) 
Jo 

where g{rj) and g(r]) = Jq 1 dxg(x; rj) are the solutions of the RG equations of the previous 
section. 

Consider first what was the traditional (replica-symmetric) situation. The RS inter- 
action parameters g(£) and g(£) are arriving to the fixed point values g* and g* (which 
are of the order of e), and then for the dependence of the renormalized mass r(£), 
according to (ft.3|) , one gets: 

r(0=r exp{A r O (3.4) 

where 

A T = 2-[(2+p)g m +pg m ] (3.5) 

At scale £ c , such that t(£ c ) is getting of the order of one, the system gets out of the 
scaling region. Since the RG parameter is by definition £ = In R, where R is the spartial 
scale, this defines the correlation length R c as a function of the reduced temperature r . 
According to ( |3.4|) , one obtaines: 

R c (r ) ~ r » (3.6) 

where v = 1/A T , eq.( |3.5| ), is the critical exponent of the correlation length. 

Actually, if the starting value of the disorder parameter g(£ = 0) = u is much smaller 
than starting value of the pure system interaction g{£ = 0) = go, the sutuation is a little 
bit more complicated. In this case the RG flow for first arrives to the pure system 
fixed point gi pure \ as if the disorder perturbation does not exist. Then, since the pure 
system fixed point is unstable with respect to the disorder perturbations, at scales bigger 
than sertain disorder dependent scale £ u the RG trajectories are eventually arriving to 
the stable (universal) disorder induced fixed point (g t , g*). According to the traditional 
theory || it is known that £ u ~ - In -. The corresponding spartial scale is R u ~ u~ u ^ a , 
and it is big it terms of the small parameter u. 

Coming back to the scaling behaviour of the mass parameter r(£), eq. (|3.4|) , we 
see that if the value of the temperature To is such that r(£) is getting of the order 
of one before the crossover scale £ u is reached, then for the scaling behaviour of the 
correlation length (as well as for other thermodynamic quatities) one finds essentially 
the pure system result R c (tq) ~ r U{pure \ However, the pure system critical behaviour 
is observed only untile R c << Ru, which imposes the restriction on the temperature 
parameter: r >> u 1 ^ = r u . In other words, at temperatures not too close to T c , 
t u << t << 1, the presence of disorder is irrelevant for the critical behaviour. 

On the other hand, if r << r u (in the close vicinity of T c ), the RG trajectories for 
<?(£) and g>(£) are arriving (after crossover) into a new (universal) disorder induced fixed 
point (g*, g*), and the scaling of the correlation length (as well as other thermodynamic 
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quantities), according to eqs.( |3.6| )-(|3.5|), is getting to be controled by a new universal 
critical exponent v which is defined by the RS fixed point (g*, g*) of the random system. 

Consider now what is the situation if the RSB scenario takes place. Again, if the 
disorder parameter u is small, in the temperature interval t u « r << 1, the critical 
behaviour is essentially controlled by the pure system fixed point, and the presence of 
disorder is irrelevant. For the same reasons as discussed above, the system gets out of 
the scaling regime (r(£) is getting of the order of one) before the disorder parameters 
start "pushing" the RG trajectories out of the pure system fixed point. 

However, at temperatures r << r u the situation is getting completely different from 
the RS case. At scales £ >> £ u (although still £ << ~ -) according to the solutions 
( |2.17| ), ( f2.22| ) the parameters g(£) and g(x;£), does not arrive to any fixed point, and 



they keep evolving as the scale £ increases. Therefore, here, according to eq. ( |3.3|) , 
the correlation length (defined, as usual, by the condition that the renormalized r(£) is 
getting of the order of one) is getting to be defined by the following non-trivial equation: 



21ni? c - r c dr ] [{2+ P )g{r ] )+pg(rj)}= In- (3.7) 

JO T 

Thus, as the temperature is getting sufficiently close to T c (in the disorder dominated 
region r << r u ) there will be no usual scaling dependence of the correlation length (as 
well as other thermodynamic quantities) like in the eq.( |3.6|) . 

Finally, as the temperature parameter r is getting smaller and smaller, what happens 
is that at scale = lni?* ~ ^ we are entering into the strong coupling regime (such 
that the parameters g(£) and g(x; £) are getting non-small), while the renormalized mass 
t(£) remains still small. 

According to the solution obtained in Appendix A, the integrals /g* drjg(r]) and 
Jo* dng(n) = G c have a finite (depending on the initial conditions) value. Thus, accord- 
ing to eq.(|3.7|), for the crossover temperature we get: 

const 

r* ~ exp( ) (3.8) 

u 



In the close vicinity of T c at r << t* we are facing the situation that at large scales the 
interaction parameters of the asymptotic (zero-mass) Hamiltonian are getting non-small, 
and the properties of the system can not be anylized in terms of simple one-loop RG 
approach. Nevertheless, the qualitative structure of the asymptotic Hamiltonian makes 
it possible to argue that in the temperature interval r << t* near T c the properties of 
the system should be essentially SG-like. The point is that it is the parameter describing 
the disorder, g{x;£), which is the most divergent. 

In a sense, here the problem is qualitetively reduced back to the original one with 
strong disorder at the critical point. It doesn't seem probable, however, that the state of 
the system will be described by non-zero true SG order parameter Q a b = ((p a (f>b) (which 
would mean real SG freezing). Otherwise there must exist finite value of r at which real 
thermodynamic phase transition into the SG phase takes place, while we observe only 
the crossover temperature r#, at which change of critical regime occurs. 
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It seems more realistic to expect that at scales ~ the RG trajectories finally arrive 
to a fixed-point characterized by non-small values of the interaction parameters and 
strong RSB. Then, the SG-like behaviour of the system near T c will be characterized by 
its highly non-trivial critical properties exhibiting strong RSB phenomena. 

3.2 Correlation Functions 

Consider the scaling properties of the spin-glass type connected correlation function: 



K(R) = ((mm) - (0(0)) (0(i?))) 2 = «0(O)0(i2)»» (3.9) 
In terms of the replica formalism one gets: 

1 - 

K(R) = lim £ K ab (R) (3.10) 

n-o n (n - 1) ^ b 



where 



K ab (R) = ((Uo)Mo)MR)UR))) (3.ii) 



In terms of the standard RG formalizm for the replica correlation function K ab (R) 
one finds: 

K ab (R) ~ (G (R)f(Z ab (R)) 2 (3.12) 

where 

G (R) = R~ iD ~ 2) (3.13) 

is the free-field correlation function, and in the one-loop approximation the scaling of 
the mass-like object Z ab (R) (with a ^ b) is defined by the RG equation: 

^1*^(0 = 2*^(0 (3.14) 

Here g aj tb(0 > is the solution of the corrsponding RG equations ( p.5|) -( ^6|) , £ = lni?, 

and Z ab (0) = 1. 

For the correlation function (|3.12|) one finds: 



K ah (R) ~ (G (i?)) 2 exp{4 f** ' d£g ab (0} (3-15) 

Jo 

Correspondingly, in the Parisi representation: g a ^ b (^) — > g(x; £) and K a ^ b (R) — > K(x; R), 
one gets: 

K(x;R) ~ (G (i?)) 2 exp{4 [** d£g(x; £)} (3.16) 

Jo 
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To realize the effects of the RSB more clearly consider again what was the situa- 
tion in the traditional RS case. Here (for p < 4) one finds that the interaction pa- 
rameter g a ^b(0 = M (0 arrives to the RS fixed point u* = e ^/^s , and according to 
eqs.( p.l5|) ,( |3.10| ) one obtains simple scaling: 



K rs (R) ~ R-^ D -V+ e (3.17) 
with the universal disorder induced critical exponent 

9 = e f~ P ^ (3.18) 
4(p-l) 

In the case of the 1-st ep RSB fixed point, eq.( |2.15| ), the situation is getting somewhat 
more complicated. Here one finds that the correlation function K(x; R) also have 1RSB 
structure: 

where (in the first order in e) 

K (R) ~ R~ 2 ( D ~ 2 ) = G 2 (R) 

(3.20) 

K X {R) ~ R-z(D-2)+e lrsb 

with non-universal critical exponent 6x rs b explicitely depending on the coordinate of the 
step xq\ 

4(4 -p) 
lb(p - 1) - px {8 +p) 
Since the critical exponent 9x rs b is positive, the leading contribution to the "observ- 



able" quantity K(R) = ((0(O)0(i?))) 2 , eq. fglOD , is given by K^R) 



K{R) ~ (1 - x^K^R) + x ^ (i?) ~ i?- 2 ( B - 2 )+^ (3.22) 

But the difference between the 1RSB the RS cases must be observed not only in 
the result that their critical exponents 9 of the correlation functions K(R) must be 
different. According to the traditional SG philosophy ||, the result that the scaling 
of the RSB correlation function K a b(R) or K(x; R) does depend on the replica indices 
(a, b) or the replica parameter x, eq. fl3.19Q , indicates that in different measurements of 
the correlation function for the same realization of the quenched disorder one is going 
to obtain different results, K (R) or Ki(R), with the probabilities defined by the value 

of Xq. 

In real experiments, however, one is dealing with the quantities averaged in space. In 
particular, for the two-point correlation functions the measurable quantity is obtained 
by integration over the two points, such that the distance R between them is fixed. Of 
course, the result obtained this way must be equivalent simply to K(R), eq. Q3.22D , found 
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by formal everaging over different realizations of disorder, and different scalings K (R) 
and K\{R) can not be observed this way. 

Nevertheless, for somewhat different scheme of the measurements the qualitative dif- 
ference with the RS situation can be observed. In spin-glasses it is generally believed that 
RSB can be interpreted as factorization of the phase space into (ultrametric) hierarchy 
of "valleys", or local minima pure states separated by macroscopic barries. Although 
in the present case the local minima configurations responcible for the RSB can not be 
separated by infinite barriers, it would be natural to interpret obtained phenomenon as 
effective factorization of the phase space into a hierarchy of valleys separated by finite 
barriers. Since the only relevant scale in the critical region is the correlation length the 
maximum energy barriers must be proportional to R®(t), and they are getting divergent 
as the critical temperature is approached. In this situation one could expect that be- 
sides the usual critical slowing down (corresponding to the relaxation inside one valley) 
qualitatively much bigger relaxation times would be required for overcoming barriers 
separated different valleys. Therefore, the traditional measuremets of the observables 
in the "thermal equilibrium" can actually correspond to the equlibration within one 
valley only and not to the true thermal equilibrium. Then in different measuremets (for 
the same sample) one could be effectively "trapped" in different valleys and thus the 
traditional spin-glass situation is restored. 

To check whether the above speculations are correct or not, like in spin-glasses one 
can invent traditional "overlap" quantities which could hopefully reveal the existance of 
the multiple valley structures. For instance, one can introduce the spartially averaged 
quantity for pairs of different realizations of the disorder: 

Kij(R) = d D r{<j>(r)<f>{r + tf)>,(0(r)0(r + R)) 3 (3.23) 

where i and j label different realizations, and it is assumed that the measurable thermal 
average corresponds to a particular valley, and not to the true thermal average. If the 
RS situation takes place (so that only one global valley exists), then for different pairs of 
realizations one will be obtaining the same result ( |3.17| ). On the other hand, in the case 
of the 1RSB, according to the general theory of the RSB ||, after obtaining statistics 
over pairs of realizations for Kij(R) one has to be getting the result K (R) with the 
probability x , and Ki(R) with the probability (1 — x ). 

Consider finally what would be the situation if a general type of the RSB takes place. 
According to the qualitative solution ( [2.17| )-( [27L8| ), the function g(x;£,) does not arrive 



to any fixed point at scales £ >> ~ -In-. Therefore, at the disorder dominanted 
scales R >> R u ~ u~ v ' a » 1 there must be no scaling behaviour of the correlation 
function K(R). Near the critical scale ~ 1/u the qualitative behaviour of the solution 
g{x;£) is shown in eq. ([2.17|) . Therefore, according to eq. ( |3.16| ), near the critical scale 



~ exp(l/w) for the correlation function K(x; R) one obtaines: 

^-2(0-2)^ _ MllljR )~4a = K X {R)- for (1 - x) « A(R) 
K(x;R)~{ (3.24) 

R -2{D-2) = Q2( R j = Kq . for (1 _ x ) » A (i?) 
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where A(R) = (1 - ulnR) — > u « 1 as R — > 

At the critical scale one has (1 — ulni?*) ~ u, and according to eq. (p.24j) the shape 
of the replica function K(x; R) must be " quasi- lstep": 

f u- 4a exp{-^^} = K{\ for (1 - x) « u 
K(x; RJ ~ J " (3.25) 

1 expj-^^} = Jf*; for (1 - x) >> u 

According to the above discussion of the observale quantities for the 1-step RSB 
case, the result ( |3.25| ) could be measured for the spartially averaged overlaps of the 
correlation functions K i j{R) 1 eq. (|3.23|) , for the statistics of pairs of realizations of the 



disorder. Then, for the correlation function Kij(R) one is expected to be obtaining the 
value K\ with the small probability u and the value K$ with the probability (1 — u). 
Although both values K* and Kq are expected to be exponentially small, their ratio 
Kl/K* ~ u~ 4a must be big. 

Finally, at scales R » R* we are entering into the strong coupling regime, where 
simple one-loop RG approach can not be used any more. 

3.3 Specific Heat 

According to the standard procedure the leading singularity of the specific heat can be 
calculated as follows: 



C ~ J d D R[mo)<j?(R)) - <0 2 (O))(0 2 (it>)>] (3.26) 
In terms of the RG scheme for the correlation function: 



W(R) = mOW(R)) - <<W)> mR)) (3-27) 

one gets: 

W{R) = (G Q (R)) 2 m 2 {R) (3.28) 

where Gq(R) = R~( D ~ 2 ^ is the free field two-point correlation function, and the mass- 
like object m(R) is given by the solution of the following (one-loop) RG equation (c.f. 

eq.Q) : 

d n 

-lnm(0 = -[(2+^(0-^^(0] (3-29) 

Here, as usual, £ = \nR, and the renormalized interaction parameters g(£) and g a ^b{€) 
are the solutions of the replica RG equations ( p.5|) -( ^6|) . In the Parisi representation, 

g^tb(0 -> g(x;0, one g ets: 

m{R) =exp{-(2+p) f^d^g^-p [ * d£ /* dxg{x;C)} (3.30) 

Jo Jo Jo 

Then, after simple transformations for the singular part of the specific heat, eq.( |3.26| ), 

one gets: 
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C~ / d£exp{e£-2(2 + p) / dr l ~g(r l )-2p dng{r,)} (3.31) 

JO JO Jo 

where ~g(rj) = Jq 1 dxg(x; rj). The infrared cut-off ^ max in ( |3 . 3 1| ) is the scale at which the 
system get out of the scaling regime. 

Usually ^max is the scale at which the renormalized mass r(£), eq. ( |3.3j ), is getting 
of the order of one, and if the traditional scaling situation takes place, one finds that 
£ max ~ ln(l/r ). 

Again, consider first what was the situation in the traditional RS case. Here at scales 
£ >> ~ ln(l/u) (which correspond to the temperature region r << r u ~ u u l a ) 
the renormalized parameters g(n) and g(£) are arriving into the universal fixed point 
9* = e i6(p-i) ' 9* = e i6(p-i) ' ( see Section 2, eq. (|2.14j) ) and according to (|3.31|) for the 



singular part of the specific heat one finds @,||]: 

/■ln(l/r ) e f-P 

C(t ) ~ / d£exp{£[e - 2(2 + p% - 2pg*]} ~ r ^ (3.32) 

So that in the close vicinity of T c one would expect to observe new universal disorder 





induced critical behaviour with negative specific heat critical exponent a = — e 4 ^_ p 1 - ) 



(unlike positive a in the corresponding pure system). 

Similary, if the scenario with the stable 1-step RSB fixed points takes place, then 
one finds that the specific heat critical exponent a(xo) is getting to be non- universal, 
explicitely depending on the coordinate of the step x 0: 

a(x ) = -l e ( 4 -^)( 4 -^°) (3.33) 
v 07 2 16(p- 1) -px (p + 8) K ' 

In the general RSB case the situation is getting completely different. Here in the 
disorder dominated region r* << t << u u ^ a (which corresponds to scales £ u << £ << 
£*) the RG trajectories of the interaction parameters and p(£) does not arrive to 



any fixed point, and according to eq.(|3.32|) one finds that the specific heat is getting 



to be a complicated function of the temperature parameter tq which does not have the 
traditional scaling form. 

Finally, in the SG-like region in the close vicinity of T c , where the interaction param- 
eters g and ~g are getting finite, one finds that the integral over £ in eq.( |3.3ip is getting 



converging (so that the upper cut-off scale £ max is getting irrelevant). Thus, in this case 
one obtaines the result that the "would be singular part" of the specific heat remains 
finite in the temperature interval ~ around T c , so that the specific heat is getting 
non- singular at the phase transition point. 



4 Marginal case p = 4 

In the systems with the number of spin components p = 4 (in which the pure system 
specific heat critical exponent a = 0) the disorder appears to be marginally irrelevant in a 
sense that it does not change the critical exponents. Nevertheless, the critical behaviour 
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(described in terms of the logarithmic singularities) is effected by the disorder, and 
moreover, the RSB phenomena appear to relevant in this case as well. 

Consider first the replica symmetric situation: g(x; £) = g(£). For the RG equations 
(EnOD , flmD one gets: 

| = (e - 12~g)g - 4g* 

(4.1) 

f = (e - 12~g)g + Ag 2 
In the pure system (g = 0) the fixed point is: 

9pure = (4-2) 

Using eq. (|3.32 ) for the singular part of the specific heat of the pure system one easily 
finds: 

C pure (r) ~ ln(-j-) (4.3) 

Thus, although the specific heat critical exponent of the pure system is zero, the specific 
heat is still divergent in the critical point. 

For the system with disorder the (replica symmetric) asymptotic solution of the 
eqs.(O) is: 



0(0 * ir 1 - 

(4.4) 

$(0*s c + 0(£) 

where 

?(£) ~ T 2 - (4.5) 

In this case the renormalized parameters are asymptotically approaching the pure system 
fixed point g = e/12, g = (so that the disorder is marginally irrelevant). Nevertheless, 
due to slow power-law approach to the fixed point the logarithmic singularity of the 
specific heat is changing into another universal type. From the general expression (|3.31| ) 
for the singular part of the specific heat one obtains: 



C 



f HlM ^exp{ p d V [e - 120(77) - 80(77)]} (4.6) 
Jo Jo 



Using the result ( |4.4|) one easily finds: 

Crs(r) - ^ (4.7) 

One can also easily check that (unlike the systems with p < 4) the crossover from the 
pure system critical behaviour, eq.( [4.3|) , to the disorder induced one, eq. Q4.71) , takes 
place in the exponentially small temperature interval near T c : 

r u ~ exp(--) (4.8) 
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Consider now the effects of the RSB. The analytic solution of the RG equations 
(|2~T0|) , (|2~TT|) (see Appendix B) shows that there is no strong coupling regime in the 
p = 4 case, and the asymptitic behaviour (at scales £ >> 1/u) of the renormalized 
parameters can be found exactly: 



1 f l-x)« 1 



m - 12 + 9(0 



(4.9) 



(4.10) 



g(0 ~ r 3/2 - 

Here 7 = g' (x = 1) ~ u is the derivative of the starting RSB function go(x) at x — 1. 

Like in the RS case the renormalized parameters are asymptotically approaching 
the pure system fixed point g = e/12, g(x) = 0. Nevertheless, the structure of the 
asymptotic solution for the renormalized function g(x;£) near this fixed point exhibits 
strong RSB. 

However, the specific heat appears to be not effected by the RSB. According to 
eq. (|3.31| ) the leading singularity of the specific heat is defined by the integral Jq 1 dxg(x; £) = 
J}(£) and not the function g(x;£) itself. It can be shown (see eq.(B.12)) that in the 
asymptotic regime the value of g(£) coinsides with the RS asymptotics 



^(0 ~ \c l (4.11) 

Therefore, for the specific heat singularity one obtains the result coinsiding with the RS 
one, eq. (|4.7p . 

On the other hand, the asymptotic behaviour of the correlation functions, appears 
to be quite different from the results of the traditional RS solution. In the RS case, 
eq.(|4.4|), according to eq.( |3.16| ) for the correlation function 



k(r) = momRW (4.12) 

one easily finds the following result: 

K(R) ~ (G ( J R)) 2 exp{4 P* d£g(Z)} = (G (R)) 2 \nR (4.13) 

Jo 

Therefore, in the RS case the disorder provides only the logarithmic correction to the 
correlation function. 

In the case of the RSB solution, eq. Ql.9D , according to eq. (|3.16|) for the replica cor- 
relation function K(x\ R) one easily finds: 



K(x-R) ~ { ( G o W ) 2 exp{ (const) v^hTR}; (1 - x)^j\nR » 1 

1 (Go(R)) 2 ; (1 -x)y/iWR« 1 l ' ; 

Correspondingly, for the "observable" correlation function, eq.( 4.12j ), one eventually 
obtains: 
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K(R) = jf 1 dxK(x; R) ~ (G (R)) 2 exp{ (const) ^7 In R} (4.15) 
This result is essentially different from the RS one, eg. ( |4.13| ). 

5 Discussion 

In this section we summarize our conclusions concerning the random bond p-component 
Heisenberg ferromagnet and discuss the remaining issues. 

Spontaneous replica symmetry breaking coming from the interaction of the fluc- 
tuations with the multiple local minima solutions of the mean-field equations has a 
dramatic effect on the renormalization group flows and on the critical properties. In 
the systems with the number of spin components p < 4 the traditional RG flows at 
dimensions D = 4 — e, which are usually considered as describing the disorder-induced 
universal critical behavior, appear to be unstable with respect to the RSB potentials as 
found in spin glasses. For a general type of the Parisi RSB structures there exists no 
stable fixed points, and the RG flows lead to the strong coupling regime at the finite 
scale R* ~ exp(l/u), where u is the small parameter describing the disorder. Unlike 
the systems with 1 < p < 4, where there exist stable fixed points having 1-step RSB 
structures, eq.( |2.15| ), in the Ising case, p — 1, there exist no stable fixed points, and any 
RSB interactions lead to the strong coupling regime. 



If there is RSB in the fourth-order potential, one could identify a phase with a 
different symmetry than the conventional paramagnetic phase, and thus there would 
have to be a temperature Trsb & t which this change in symmetry occurs. Actually, the 
RSB situation is the property of the statistics of the saddle-point solutions only, and it 
is clear that for large enough r there must be no RSB. Therefore, one can try to solve 
the problem of summing over saddle-poin solutions for arbitrary r, aiming to find finite 
value of t c at which the RSB solution for this problem disappears. 

Of course, in general this problem is very difficult to solve, but one can easily obtain 
an estimate for the value of r c (assuming that at r = the RSB situation takes place). 
According to the qualitative study of this problem in the paper || , the RSB solution can 
occur only when the effective interactions between the "islands", (where the system is 
effectively below T c ) are getting non-small. The islands are the regions where 5r(r) > r. 
According to the Gaussian distribution for Sr(r), the average distance between them 
must be of the order of exp[— t 2 /u], so that the islands are getting distant at r > yfu. 
The interaction between the islands is exponentially small in their separation. Therefore 
at t > y/u they must be getting weakly interacting, and there must be no RSB. 

Note now that the shift of T c with respect to the corresponding pure system is also 
of the order of y/u. On the other hand, the existence of local solutions to the mean- 
field equations remines the Griffith phase []I2] which is claimed to be observed in the 
temperature interval between T c of the disordered system and T c of the corresponding 
pure system. On these grounds it is tempting to associate the (hypothetical) RSB 
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transition in the statistics of the saddle-point solutions with the Griffith transition. 
Correspondingly, it would also be natural to suggest that discovered RSB phenomena in 
the scaling properties of weakly disordered systems could be associated with the Griffith 
effects 



The other key question which remains unanswered, is whether or not the obtained 
strong coupling phenomena in the RG flows could be interpreted as the onset of a kind 
of the spin-glass phase near T c . Since it is the RSB interaction parameter describing 
disorder, g(x;£), which is the most divergent, it is tempting to argue that in the tem- 
perature interval r << ~ exp(— 1/u) near T c the properties of the system should be 
essentially SG-like. 

It should be stressed, however, that in the present study we observe only the crossover 
temperature r*, at which the change of the critical regime occurs, and it is hardly possible 
to associate this temperature with any kind of phase transition. Therefore, if the RSB 
effects could indeed provide any kind of true thermodynamic order parameter, then this 
must be true in a whole temperature interval where the RSB potentials exist. 

The true spin-glass order (in the traditional sense) arises from the onset of non-zero 
order parameter Q a b{x) =< (j) a (x)(j)b(x) >;a b, and, at least for the infinite-range 
model, Q a b develops the hierarchical dependence on replica indices obtained by Parisi 
[ |T6|j . In the present problem we only find that the coupling matrix g a b for the fluctuating 
fields develops strong RSB structure and its elements are getting non-small at the finite 
scale. Therefore, it seems more realistic to interprete discovered RSB strong coupling 
phenomena in the RG just as a new type of the critical behaviour characterized by 
strong SG-effects in the scaling properties rather then in the ground state. 

In spin-glasses it is generally believed that RSB phenomenon can be interpreted 
as a factorization of the phase space into (ultrametric) hierarchy of "valleys", or local 
minima pure states, separated by macroscopic (infinite) barries |J. Although in the 
systems considered here the local minima configurations responcible for the RSB are 
not likely to be separated by infinite barriers (otherwise it would mean true SG freez- 
ing), it would be natural to interpret obtained phenomenon as effective factorization 
of the phase space into a hierarchy of valleys separated by finite barriers. Since the 
only relevant scale in the critical region is the correlation length the maximum energy 
barriers must be proportional to R®(t), and they are getting divergent as the critical 
temperature is approached. In this situation one could expect that besides the usual 
critical slowing down (corresponding to the relaxation inside one valley) qualitatively 
much bigger (exponentially large) relaxation times would be required for overcoming 
barriers separating different valleys. Therefore, the traditional measuremets (made at 
finite equilibration times) can actually correspond to the equlibration within one valley 
only, and not to the true thermal equilibrium. Then in a close vicinity of the critical 
point different measuremets of the critical properties of e.g. spatial correlation func- 
tions (in the same sample) would exhibit different results as if the state of the system 
is getting effectively "trapped" in different valleys, and thus the traditional spin-glass 
situation will be observed. 
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Appendix A: The asymptotic solution for p < A case 

In this Appendix we derive the asymptotic solution of the RG eqs.( p.l0[) , (|2.11|) : 

j-g(x) = (e-(A + 2p)g)g(x)+Ag 2 (x)-2pg(x) dyg(y)-p^ dy(g(x) - g(y)) 2 (Al) 

—g = eg - (8 + p)g 2 + pip (A2) 

(where g 2 = J dxg 2 (x)) for the number of components p < A. 

It can be shown aposteriori that the term (e — (4 + 2p)g)g(x) in the eq.(Al) is 
irrelevant in the asymptotic regime. So, consider the equation: 

—g(x) = Ag 2 (x) - 2pg(x) / dyg(y) -p dy(g(x) - g{y)f (A3) 
a£ Jo Jo 

After taking derivative over x and after simple transformations one gets: 

^g'(x) = 2pg'(x)[(X - l)g(x) - £ dy(l - y)g'(y)] (AA) 
where A = A/p > 1. Let us introduce: 

V(x)= [ 1 dy(l-y)g'(y) (Ah) 



(A6) 



According to this definition one has: 

9\x) = -j^V'(x) 



g( x ) = J*dyg>(y) = -J*dy^-V>(y) 



Here for simplicity we consider the case g(x = 0) = (the behaviour of the solution for 
g(x = 0) 7^ in the asymptotic regime can be shown to be qualitatively the same). 
Then, for the eq.([5|) after simple transformations we get: 

^-V'(x) = -2pV'(x)[ f X dy^-V'iy) +g(0) (At) 
at, Jo 1 — y 
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where g(£) = /„* dxg(x, f ) = Jq dx(l - x)g'(x) = V(x = 0; f )• 
Let us define now: 



W(x;0= f X dy^—^V'(y) (A8) 
Jo l-y 



or 

V'(x) = \^-W'(x) (A9) 



From eq.(A7) one gets: 



-^W'(x) = -2pW'(x) [W(x) + g(0] (A10) 



Integrating over x yields: 

^-W{x) = - P W 2 (x) - 2pW(x)g(£) (All) 
at, 

(Here the integration constant is zero because W(x = 0) = 0). This equation can be 
easily solved for any given function 7?(£): 



Wpjx) exp[-2p J* dr]g(rj)} 

1 + pW (x) fl> dt exp[-2p fl drjg(rj)} 

where: x 

W o (x) = W(x;Z = 0) = - f dy{\-y)g' Q {y) (A13) 



and go(x) = g(x;t; = 0). Coming back through the definitions (A8) and (A5) for the 
function g(x;^) one gets: 



o * [i-p js d v e( v ) jg dz(\ - z) g ' {z)Y 



g(x;t)= / dy , M^Z^ {AU ) 



where: 

9(0 = exp[-2p f dr,g(r,)] (A15) 
Jo 

Integrating f^ dxg(x; £) = one gets the equation for the unknown function #(£): 

Now the problem is to find the asymptotic behavior of <?(£). 
Let us introduce: 

G(0= (A17) 

Integrating (A16) we obtain: 

Jo y [l-M(0/ w dz(A-zy o (*)] 
24 



where: 



Let us redefine 



ft 

A(£) = / dr)ex.p[-2 P G(r))} (A19) 
Jo 



m = (A(o)- 1 = -j- — - - - (A20) 

Then: 

Now, let us redefine again: 

m=P [ 1 dy(X-y)g > (y) + m (A22) 
Jo 

From eq.(A21) we get: 

A) [p^ dz(A - 3)^(2) +0(0] 
Assuming that 0(0 is small, (A23) can be estimated as follows: 

C(0 =G C + f dy(l-y)g' (y)[- T —- 1 — — ] (A24) 

Jo pj y dz(X- z)g' (z) + 0(0 pj y dz( y \-z)g' (z) 



or 



G(0 = G c - 0(0 [ dy ri r TT , , n VTT^ ( A25 ) 

Jo [pjy dz(X - z)g' (z) + <j){0][PJy dz(X - z)g' (z)\ 

where 

G c = f 1 dy (A26) 

For 0(0 << 1 the leading contribution in the integral in (A25) comes from the 
vicinity of y = 1. Assuming that g' (y = 1) = 7^0, this contribution can be estimated 
as follows: 



G(0 * G c - 0(0 X 1 . 



(i-y)7 



[P7(A-1)(1-!/)+<K£)]p7(A-1)(1-2/) ~~ 

(A27) 



Such that, as — > 0, the value of G(0 goes to finite value G c , but near this point the 
behavior of this function is non-analytic. 



25 



Now let us assume that there exists sertain scale £ c , such that 0(£ — ■> £ c ) — * 0, and 
consider the behavior near £ c . Coming back to the definition (A20) we can estimate: 

m = ltidr ] exp(-2pG(r ] ))]-t = 
= dr)eM-2pG(v)) ~ ft dr ] exp(-2pG(r ] )))- 1 ~ 
~ (ft d V exp(-2pG( V )) - exp(-2pG c )^ c - ^ 



(A28) 



exp(-2p Ge ) (gc _ g) 



J ( f c dr,cxp(-2pG(»,)) [/^ ^oxp(-2 P G(r,))] 2 

Comparing this result with (A22), we find that: 

0(0 ~ a(£ c - (^29) 
where the parameters £ c and a are defined by: 

1 f 1 



f£ c dr}exp(-2pG(r})) 

and 



pf dy(X-y)g' (y) (A30) 
Jo 



a = ex P( 2 ^ = \p f 1 dy{ X - y)g' {y)f exp(-2 P G c ) (A31) 
[f^driexp(-2pG(r]))] 2 Jo 

Let us estimate the parameters £ c and a by the order of magnitude. The characteristic 
value of the initial function go(x) is of the order of u « 1, which is the characteristic 
value of the quenched disorder. If the initial function go(x) does not have special anomaly 
near x — 1, then its derivative 7 must also be of the order of u. Then, the above integrals 
can be estimated as follows: 

Jo p f y dz(X - z)g' (z) 

[ 1 dy(\-y)g' (y)^u (A33) 
Jo 

drjexp(-2pG(rj)) ~ £ c 0434) 

Jo 

Thus, from (A30) and (A31) for the parameters £ c and a we find: 

£ c ~ - (^435) 

■u 

a ~ w 2 (A36) 

Now we can describe the qualitative behavior of the asymptotic solution. According 
to (A27): 
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(A37) 

Therefore the value of the intergal Jq 1 dxg(x; £) = g(£) is formally getting divergent at 
finite scale £ c ~ l/u. 

Coming back to the result (A14) for the function g(x; £) we have: 

g ( x; f ) = 9(0 f* dy ^ 

= [ d 1 1 f x dv 



(A38) 



Therefore, when approaching the critical scalej £ — > £ c , the values of g(x; £) are formally 
getting big in the narrow interval (1 — x) « A(£), where: 



A(O~^(&-0~(1-«O (^39) 



In this interval: 



- ° X 1 . d V [p(A-l) 7 (l4)+a(«c-?)] 2 



(A40) 
i i 



p(A-l) £ c -£ 



where a ~ 1 is a (non- universal) constant. 

Therefore, the considered RG approach can be applied only up to the scales, such 
that (1 — u£) ~ u (until the value of the parameter g\ is getting non small). 

Now let's come back to the equation for the diagonal parameter g (A2). According 
to the asymptotics obtained above we can estimate the value of g 2 . Since the leading 
contibution comes from the region A(£) near x — 1, we get: 

? =/; <fa9 '(*i{)~(i-«£)( i ^5 i =«* T ^ (^4D 

Therefore, from the equation (A2) we see that g diverges as the logarithm: 
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S(£)~ U ]n— L- (AA2) 
1 — ut, 

Thus, in the region near x — 1 where the value of g(x\ £) was obtained to be of the order 
of uj (1 — «£) the first term (e — (4 + 2p)g)g(x) in the eq.(Al) is much smaller than the 
other terms: 

1 u 2 

99(x) ~ ~, 7 In z t « 71 ~ <? 2 (x) (A43) 

Of course, the above asymptitic solution does not make possible to obtain the be- 
havior of the function g(x;£) in the whole interval [0,1] for all £. Nevertheless, the 
numerical solution of the general RG equations (Al), (A2) clearly demonstrates that at 
large scales the function g(x;£) qiuckly goes to zero for all x not too close to 1, while 
in the narrow region near x = 1 the values of this function are getting divergent. Thus, 
the behaviour of the asymptotic solution for g(x; £) in the vicinity of the critical scale 
£ c could be qualitatively represented as follows: 

f a^; for (1 - x) « A(£) 
g(x;£)~l (AM) 
{ 0; for (1 - x) » A(0 

where A(£) = (1 — «£) —>■ u << 1 as £ — > £ c , and a is a positive non-universal constant. 

The obtained asymptotics can also be easily generalized for the situation when g(x = 
0) 7^ 0. One has to write: g(x;£) = (the obtained solution) +g(x = 0;£), then put it 
into the equation, obtain the equation for g(x = 0;£), and find the asymptotics for 
g(x = 0;£) It's straightforward to check that qualitatively it dosn't change the above 
results. 



Appendix B: The asymptotic solution for p = 4 



In p = 4 case the asympotic solution of the equations ( |2. 10| )-( pTTT| ) can be obtained 
as follows. 

Redefining the diagonal parameter g(£): 



we get: 



^-g{x) = -12q(0g(x) + Ag 2 (x) - 8g(x) C dyg(y) - 4 f dy(g(x) - g{y)) 2 (B2) 
at; Jo Jo 

^(0 = -eg-12g 2 + 4^ (B3) 
Then, proceeding like in the Appendix A, instead of eq.(A7) we obtain: 
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±V'(x; = SV'(x; £)V(x; - 12y'(x; 0?(0 (54) 
Integration over x yields: 

±V(x; = -4V 2 (x; - 12q(t)V(x; (55) 

(the integration constant is zero, since V(x = 1) = 0). The solution of this equation for 
any given function q(£) is: 



v( X] o = ^\ 7; r r,; ^ (56) 



Vft(a;) exp{-12 J| driqfy)} 
1 + 4\/ (x) j* ^ exp{-12 cftg(*)} 

where V (x) = V(x; £ = 0) = £ - y)g' (x). 
Coming back to the function g(x; £) we get: 



g(x; = f <%'(</) + y(x = 0; = - f dy-^V'(y) + g(x = 0; (57) 
Jo jo 1 — y 

Using (B6) we find: 

« = y * TTZIf^ "til r> FTvui + ^ = 0: « < B8 > 

Jo [l + 4/ 2/ i ^(l-z)^(z)/ ^77exp{-12/ ?? ^(i)}] 2 

Putting this result back into the original equation (B2) we get the equation for g(x — 
0:0: 



— <7(x = 0; = -12q{S)g{x = 0; - 4^ 2 (a; = 0; - Sg(x = 0; 0^(0 (59) 
a? 

where 

a(£) = f 1 dx T dv ^ (t/)exp{-12/ € ^g(77)} 

io 7o y [l + 4/ s/ 1 ^(l-^ (^)/ ^r ? exp{-12/J'^(t)}] 2 

Let us assume now that the parameter q(£) decays as ~ with s > 1. Then the 
integral P drjq(r]) is converging at large £, and for the exponent in (B8) we find that it 
is equal to a constant of the order of one: exp{— 12 ffj drjq(rj)} = A. 

Correspondingly, instead of eqs.(B8),(B10) we get: 

[dy j— ^jrwrir \ , 12 +^ = Q;0 (5ii) 

Jo [1 + 4^/j, dz(l -z)g' (z)] 2 

and 
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g(0 ~ f} dx /* dy 

(512) 

1+4A?3 

where # = Jq dxg (x) = ft dx(l - x)g' (x). 

Simple analycis of the integral in eq.(Bll) shows that actually it is the non-zero 
derivative g' (x) near the point x = 1 which is important in the asymptotic regime. 
Whatever the function g (x) is in the region (1 — x) » (7O ^ is always decaying 
like £~ 2 there, while for (1 — x) « (7£)~ 1/2 the decay is (7£)~ 1/2 > where 7 = g' (x = 1): 



r 2 , (i-x)» 1 



. n ^ _ 'i 7 ' (513) 



— ^ (1 - x ) « —r= 

Besides, using (B12), from eq.(B9) one finds that y(x = 0;£) ~ £~ 2 . 

Note now that according to the above asymptotic behaviour of the function g(x;£) 
at scales £ >> \ ju the leading contribution to the quantity g 2 = Jq 1 <? 2 (:c; £) comes from 
the region (1 — x) « 



7 ~ 7t ( 7t r = r3 " (BU) 



Then, coming back to the eq.(B3) we find: 



q(0 ~ exp(-eO J dtr 3/2 exp(+et) ~ ai r 3/2 + a 2 C 5/2 + - ~ T 3/2 (515) 

which is selfconsistent with the assumption g(£) ~ (with s > 1) made above. 

Therefore, the asymptotic behaviour of the solution for p = 4 at scales £ >> 1/tt is 
given by eq.(B13). 
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